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Abstract 

We calculate the spectral gap of the Markov matrix of the 
totally asymmetric simple exclusion process (TASEP) on a ring 
of L sites with N particles. Our derivation is simple and self- 
contained and extends a previous calculation that was valid only 
for half-filling. We use a special property of the Bethe equations 
for TASEP to reformulate them as a one-body problem. Our 
method is closely related to the one used to derive exact large 
deviation functions of the TASEP. 
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1 Introduction 

The asymmetric simple exclusion process (ASEP) is a driven diffusive system 
of particles on a lattice interacting through hard-core exclusion, that serves 
as a basic model in various fields ranging from protein synthesis to traffic flow 
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(for a recent review, see Schiitz 2001). In non-equilibrium statistical physics, 
the ASEP plays the role of a paradigm thanks to the variety of phenomeno- 
logical behavior it displays and to the number of exact results it has led to 
in the last decade (see, e.g., Derrida 1998). In particular, ASEP is an inte- 
grable model, i.e., the Markov matrix that encodes its stochastic dynamics 
can be diagonalized by Bethe Ansatz, as first noticed by Dhar (1987). Thus, 
the spectral gap of the Markov matrix, i.e., the difference between the two 
eigenvalues with largest real parts, that characterizes the longest relaxation 
time of the system, can be calculated exactly: this was first done for the 
totally asymmetric simple exclusion process (TASEP) at half filling (Gwa 
and Spohn 1992) and later, using a mapping into the six vertex model, Kim 
(1995) treated the case of the general ASEP at arbitrary filling. In both 
works, the calculations are complicated though the final result for the gap is 
fairly simple. In a recent work (Golinelli and Mallick 2004), we presented a 
concise derivation of the TASEP gap at half filling that circumvents most of 
the technical difficulties thanks to an analytic continuation formula. How- 
ever, the half filling condition seemed to play a crucial role in our derivation 
(as well as in the calculation of Gwa and Spohn). Our aim in the present 
work is to show that our method can be extended to the arbitrary filling 
case. 

We shall study the TASEP on a periodic one-dimensional lattice with L 
sites (sites i and L + i are identical). The TASEP is a discrete lattice gas 
model on which each lattice site % (1 < i < L) is either empty or occupied 
by one particle (exclusion rule). Particles evolve according to stochastic 
dynamical rules: a particle on a site % at time t jumps, in the interval between 
times t and t + dt, with probability dt to the neighbouring site % + 1 if this 
site is empty. As the system is periodic, the total number N of particles is 
conserved and the density (or filling) is given by p = N/L. A configuration 
of the TASEP can be characterized by the positions of the N particles on the 
ring, (x\, X2, ■ ■ ■ , xn) with 1 < x\ < X2 < ■ ■ ■ < xn < L. If ipt(xi, ■ ■ ■ , xn) 
represents the probability of this configuration at time t, the evolution of 
tpt is given by the master equation dip t /dt = Mip t , where M is the Markov 
matrix. A right eigenvector ip is associated with the eigenvalue E of M if 
Mvfj = Eip. Thanks to the Perron- Frobenius theorem, we know that the zero- 
eigenvalue of M, which corresponds to the stationary state, is non-degenerate 
and that all the other eigenvalues of M have a strictly negative real part. In 
the stationary state, all configurations have the same probability, given by 
N\(L-N)\/L\. 
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In the next section, we present the Bethe Ansatz equations and restate 
them as a single self-consistency equation. We then calculate the TASEP 
spectral g ctp clS db function of the density p in the limit of a large system size, 
L — > oo (Section 3). 

2 The Bethe Ansatz Equations 

The Bethe Ansatz assumes that the eigenvectors ijj oi M can be written in 
the form 

1p(xx, ...,X N )= J2 A ° C(l) 4(2) • • • Z l{N) , (!) 

where is the group of the N\ permutations of N indexes. The coefficients 
{A a } and the wave-numbers {zi, . . . , zn} are complex numbers determined 
by the Bethe equations. In terms of the fugacity variables Zi = 2/^ — 1, 
these equations become (Gwa and Spohn 1992) 

N 7, - 1 

(l-^(l + Z 4 )^ = -2 L n^-T with i = l,-,N. (2) 

We note that the right-hand side of these equations is independent of the 
index i: this property is true only for the totally asymmetric exclusion process 
and not for the partially asymmetric exclusion process where the particles 
can also jump backwards. Introducing an auxiliary complex variable Y, 
the Bethe equations (j2J can be reformulated as explained below (for more 
details, see Gwa and Spohn 1992; Golinelli and Mallick 2004). Consider the 
one variable polynomial equation of degree L, 

(1- Zf (1 + Z) L ~ N = Y , (3) 

and call (Zi, Z 2 , . . . , Z{) the L roots of this equation. For a given value of Y, 
the complex numbers (Zi, Z 2 , . . . , Z L ) belong to a generalized Cassini oval 
defined by the equation 

\Z - l\ p \Z + l\ 1 ~ p = r with r=\Y\ 1/L , (4) 

p being the density of the system. The topology of the Cassini oval depends 
on the value of r (see the figure). Defining 



r c = 2 f f(l-pY 1 -"\ 



(5) 



O. Golmelli, K. Mallick — TASEP spectral gap 



-N+1 



-N 



-N-1 







7 




<l-2nX 

\^ / i f T 








,1-1 J 






J 





Figure 1: Roots of the equation (1 — Z) N (1 + Z) = Y . Here N = 6 and 
L = 15 (the filling is p = 2/5) with Y = e l ^r L for <fi = tt/2; curves are drawn 
forr/r c = 0.8 (blue), 1 (black) and 1.2 (red). See text for more explanations. 



we find that for r < r c , the locus of the Z's consists of two disjoint ovals, 
with the roots [Z\, Z 2 , . . . , Z^) belonging to the oval on the right and the 
roots (Zn + i, Zn +2 , . . . , Zl) to the oval on the left. For r = r c , the Cassini 
oval is a deformed lemniscate of Bernoulli with a double point at Z c = 1 — 2p. 
For r > r c , the oval is made of a single loop. The labelling of the L roots 
(Zi, . . . , Z L ) is shown in the figure. 

We now introduce a monotonous function c : {1, . . 
that selects N fugacities among the L roots (Zi, . . . , Zl) 
tion A C (Y) as 



,N} -* {1,...,L} 
Defining the func- 



N 



a c (y) = -2 l n 



z <j) - 1 
+ 1 



(6) 



j=l ^c(j) 

the Bethe equations become equivalent to the self- consistency equation 

A C (Y) = Y. (7) 

Given the choice function c and a solution Y of this equation, the Z^'s are 
determined from Eq. (JHJ) and the corresponding eigenvalue E c is given by 



2E r 



-N 



N 



(8) 
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The choice function Co(j') = j that selects the N fugacities Zi with the 
largest real parts provides the ground state of the Markov matrix. The 
associated A-function and eigenvalue are given by 

A (Y) = -2^n^-T> (9) 

j=l Z i + 1 

N 

2Eo = -N + Y,Zj- (10) 

3=1 

The equation A (Y) = Y has the unique solution Y = that yields Zj = 1 
for j < N and provides the ground state with eigenvalue 0. We shall also 
need the following formula in the sequel: 

To derive this identity, we raise equation (jOJ) to the iV-th power, use equa- 
tion @ and take the logarithm of the result; however an additional term, 
which is a discrete constant (a multiple of 2m /N) appears in the calculations. 
This constant vanishes identically because both terms in equation (|TT|) be- 
come real numbers in the limit Y —>■ 0. 

The spectral gap, given by the first excited eigenvalue, corresponds to the 
choice Ci(j) = j for j — 1, . . . , N — 1 and ci(N) = N + 1 (Gwa and Spohn 
1992). The associated self-consistency function and eigenvalue are given by 
(using equations ISHTUjO 

MY) = MY)^1^±1, (12) 

^N+l + 1 — 1 

2E 1 = 2E + {Z N+1 - Z N ). (13) 

Thus, from equation (|T2"j) . the self-consistency equation that determines the 
gap reads 

= ln^ = ln M0_ h fi^ . (14) 

The excitation corresponding to the choice function c(j) = j + 1 for 
j = 1, . . . , N — 1 and c(iV) = L, leads to the complex-conjugate eigenvalue 
-Ei. The eigenvalue E\ corresponding to the first excited state is obtained 
by solving the self-consistency equation Ai(Y) = Y where Ax is defined in 
equation (|12|l . We shall solve this equation in the limit N,L —> oo, keeping 
p = N/L constant. 
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3 Calculation of the first excited state 

As in (Golinelli and Mallick 2004), we start with the Taylor expansions, of 
Aq(Y) and Eq(Y) i n the vicinity of Y = and valid for arbitrary values 
of N and L. In our previous work, we considered only the half-filling case, 
L = 2N, for which the polynomial equation Q reduces to (1 — Z 2 ) N = Y 
and can be solved explicitly in terms of Y yielding the Taylor series of A and 
E . An explicit solution of equation Q can not, however, be obtained for 
an arbitrary density p. This major technical difficulty can be circumvented 
thanks to a contour integral representation similar to that used for calculating 
large deviation functions (Derrida and Lebowitz 1998; Derrida and Appert 
1999; Derrida and Evans 1999). 

When Y — > 0, the N roots (Zi, . . . , Z^) of equation Q with the largest 
real parts converge to +1, whereas the L — N remaining roots converge to 
-1. We now consider a positively oriented contour 7 that encircles +1 such 
that for sufficiently small values of Y the roots (Zi, . . . , Z^) lie inside 7 and 
(Z N+ i, . . . , Z L ) are outside 7. Let h(Z) be a function that is analytic in a 
domain containing the contour 7. We also define 



Because, by definition, the Zj are the zeroes of P(Z) = Y, we obtain, from 
the residues theorem, 



Expanding the denominator in the contour integral for small values of Y 
(thanks to the formula (P - Y)- 1 = ET=o Yk / pk+1 valid for \ Y \ < \ p \)i we 
obtain the Taylor expansion 



P(Z) = (1-Z) N (1 + Z) 



L-N 



(15) 




(16) 




N 



(17) 
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The second equality is derived by integrating by parts the terms with k > 1; 
the term Nh(l) in the third equality is obtained from the residues theorem. 
We shall need the following identity, valid for any positive integers P and Q: 

1 i 1 .dz=- 2 -Wt g r 2 V as) 



2m J 7 (1 - Z) p (l +Z)Q \ P - 1 

Using equation (fTTj) with h(Z) = In f"Mp), we obtain from equations (JTTj) 
and (HHJ) 

= £U*J^" (19) 

Similarly using equation ((T7j) with /i(Z) = Z — 1, we obtain from equa- 
tions (|nS|) and ((TBI) 

In the limit L —>■ oo and with p fixed, we obtain from the Stirling formula 



MY) 1 ^ 

y ^ J27T P (1 - p)L 



, = ■ Mfr • (2D 



where r c was defined in equation © and the polylogarithm function of index 
s, Li s , is given by 

U s (z) = —- / — = Y — . 22 

v ; r(s) Jo e*-z ^h s K J 

The function Li s is defined by the first equality on the whole complex plane 
with a branch cut along the real semi-axis [1, +oo); the second equality is 
valid only for \z\ < 1. 

The limit found in equation (f2*Tj) suggests that Y can be parameterized 
as follows: 

Y = -r^e U7T , (23) 

u being a complex number with — 1 < Im(it) < 1 and which remains finite in 
the limit L — > oo. Thus |F| 1/,L ~ r c and roots of the type Zk, ZL-k and Z^±k 
where k is a fixed positive integer are close to the lemniscate double-point 
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Z c = 1 — 2p. Writing Z = 1 — 2p + 2£, with (<1, and taking the logarithm 
of equation (J2J) we obtain 

For fixed k, the values < lead to and Zl_£ , and k > 1 lead to Zj. 

and Z^+k- Neglecting terms of order (9 ^L~ 3 / 2 ), we have 



;„ = , _ 2p + Ti^SSfA (u - + - + ■ ■ ■ (25) 



Z ff+1 = 1 - 2„ + »^I*Iz>> (u + i)" 2 + f + i) + . . . (26) 

Substituting equations pijl. (J23j) and (|2Tj) in equation (JHj), we obtain the 
large L limit of the Bethe equations for the gap which, at the leading order, 
reads as 



Li 3/2 (-e^) = 2tir 



[u + i^-iu-i) 1 ' 2 



(27) 

This equation is the same as that obtained in (Golinelli and Mallick 2004) 
and its solution is given by: 

?i = 1.119 068 802 804474... (28) 

We can now calculate the eigenvalue corresponding to the first excited state. 
^From equations (fT3j) . (|T4*j) and (|2*UJl . we obtain 

2E 1 = (2E + Z N+1 - Z N ) + 2p(l - p) In ^2 = ( 29 ) 

^ P{1 ~ P) {kN) Yk 7 , oh m (1-Z N )(1 + Z N+1 ) 

^7T7 nofc/-i +(^iV+i-Z iV )-2p(l-p) In -. 

~ k(kL - l)2 kL 1 (1 + Z N )(1- Z N+1 ) 

The large L limit of this expression is found by using the Stirling for- 
mula, the expansions ()25|) and ([26)1. and the parameterization of K given by 
equation (J2*3j) . We thus obtain 



L3/2\/ 7T 

4«vr . . . 

+— (l-2p), (30) 
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where u is the solution of the equation (|27j) . its numerical value being given 
in equation (|28|) . Comparing with the value Ei p= i/2 obtained at half-filling 
we conclude that 

E x = 2y/p(l-p) E 1>p=1/2 + ^(1 - 2p) . (31) 

This equation agrees with the one derived by Kim (1995). We notice that this 
eigenvalue has a non- vanishing imaginary part when the density is different 
from one half. 

4 Conclusion 

The asymmetric exclusion process can be mapped into the six vertex model 
and is thus an integrable model that can be solved by Bethe Ansatz. Kim 
(1995) has used this technique to calculate gaps and crossover functions for 
the generic asymmetric exclusion process but the calculations are very com- 
plicated. However, for the totally asymmetric exclusion process the Bethe 
equations can be reduced to a single polynomial equation and their analysis 
becomes much simpler as shown in the present work. We have calculated 
the gap of the TASEP for arbitrary filling. Our method is closely related 
to that used to calculate large deviation functions (see Derrida 1998, for a 
review) and leads to derivations that are far simpler than the ones presented 
in previous works. This technique can be generalized to calculate any finite 
excitation close to the ground state of the system. 
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